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It would be hard to find a cosmologist today who 
does not believe that the vast bulk of the Universe 
(ninety-five percent or more) is hidden from our eyes. 
We review the evidence for this remarkable consensus, 
and for the latest proposal, that the mysterious dark 
matter consists of as many as four separate ingredients: 
baryons, massive neutrinos, new "exotic" dark matter 
particles, and vacuum energy, also known as the 
cosmological constant (A). Of these, only baryons fit 
within standard theoretical physics; the others, if their 
existence is confirmed, will mean rewriting textbooks. 
Fresh experimental evidence has recently appeared for 
and against all four components, so that the subject is 
in a state of turmoil and excitement. The past few years 
in particular have seen the fourth (vacuum) component 
come into new prominence, largely at the expense of 
the third (exotic dark matter). We conclude our review 
by exploring the possibility that the energy density of 
the vacuum is in fact so dominant as to leave little 
room for significant amounts of exotic dark matter. 
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1 The Four Elements of 
Modern Cosmology 

Hear first the four roots of all things: bright Zeus 
(fire), life-giving Hera (air), Aidoneus (earth) and 
Nestis (water), who moistens the springs of mor- 
tals with her tears. 

- Empedokles, Fragments, c. 450 B.C. 

Like the philosophers of antiquity, modern cos- 
mologists have divided the physical world into four 
different realms, each characterized by its own 
length scale and dominated by an increasingly rar- 
efied species of invisible matter which, while not 
seen directly, is inferred to exist from its gravita- 




Fig. 1. The "four elements" of modern cosmology (adapted 
from a figure in a 1519 edition of Aristotle's Libri de Caelo) 



tional influence on visible matter as well as the 
geometry of the Universe. That which is seen - 
mostly luminous hot gas, stars and galaxies — has 
a density (denoted ^lum) hundreds of times smaller 
than the estimated total density (f2 T oT) of matter 
and energy in the Universe. Complementing it is 
a considerably larger amount of so-called baryonic 
dark matter (of density B dm)) which has a simi- 
lar composition but does not shine. Even together, 
however, these two ingredients (with combined den- 
sity SI bar = ^lum + ^bdm) make up less than five 
percent of the total density. Baryonic ("ordinary") 
matter is thus relegated to a minor role in the cos- 
mic scheme, a development which has rightly been 
seen as a "second Copernican revolution," and one 
which lends a double meaning to the identification 
of baryons with "earth," the first element of the 
new cosmology (Fig. 1). 

Observations of clumps of baryonic matter on 
scales larger than the solar system (galaxies, clus- 
ters of galaxies) have led many to postulate the exis- 
tence of an additional form of matter: the "exotic" 
cold dark matter (CDM) particle, whose collective 
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density f2 CDM would far exceed that of the baryons. 
Theoretical particle physics provides several plausi- 
ble candidates, and it could even be that this com- 
ponent of the cosmic fluid — identified with "water" 
in Fig. 1 — itself has several ingredients. However, 
despite intensive searches, none of these candidate 
particles has yet been observed. 

On still larger scales — those relevant to the for- 
mation of galaxy clusters in the early universe - 
a different form of dark matter may be operative: 
the neutrino. That this particle exists is unques- 
tioned; but the extent of its contribution (Q v ) to 
the density of the Universe is not yet clear. If the 
neutrino has a small (or zero) rest mass, then it is 
always relativistic and can be treated for dynamical 
purposes like the photon. (We have therefore clas- 
sified it together with light as the "air" of the new 
world view, Fig. 1.) In this case, neutrino contribu- 
tions can be combined with those of photons (0 7 ) 
to give the total radiation density of the Universe, 
£l R = £l u + fL,. This is known to be insignificant 
at present. If, on the other hand, neutrinos are 
sufficiently massive, then they are no longer rela- 
tivistic on average, and belong instead (with bary- 
onic and exotic dark matter) under the category of 
dust-like (zero-pressure) matter, with total matter 
density ft M = B ar + ^cdm + &v Only in the lat- 
ter case do neutrinos play a significant dynamical 
role in the present Universe. Recent experimental 
evidence has been taken by many to support this 
second scenario, implying (in most models) a col- 
lective neutrino density well below that of baryonic 
matter, but (in others) a density possibly rivalling 
that attributed to exotic matter. 

Influential only over the largest of scales — the 
cosmological horizon — is the outermost species of 
invisible matter: the vacuum energy (also known by 
such names as dark energy, quintessence, x-matter, 
the zero-point field, and the cosmological constant 
A). This profusion of nomenclature betrays the fact 
that there is at present no consensus as to where 
vacuum energy originates, or how to calculate its 
energy density from first principles. Existing 

theoretical estimates of this latter quantity range 
over some 122 orders of magnitude, prompting most 
cosmologists until very recently to disregard it al- 
together. New observations of distant supernovae, 



however, suggest that the vacuum not only gravi- 
tates, but that its effective density exceeds that of 
all the other forms of matter put together. Thus 
the "standard model" of cosmology has evolved (in 
the space of three years) from one in which Q A = 
in practice (and is as large as ~ 10 122 in theory), 
to one in which there is widespread agreement that 
S1 A must be of order unity — indicating to us that 
cosmology, far from being a "solved problem" as 
reported by some authors, remains in a far-from- 
settled state. We have identified vacuum energy 
with "fire" in Fig. 1. 

Four kinds of invisible matter: dark baryons, 
exotic particles, neutrinos and vacuum energy - 
three of which imply new physics, and any one of 
which (with the possible exception of the neutrino) 
outweighs the entirety of the visible Universe. Are 
they all necessary? Our purpose in this review will 
be to re-examine the steps leading up to this con- 
clusion. We begin with the relevant cosmological 
equations in §|2|. The "four elements" are then in- 
troduced in turn: dark baryons (§||), exotic matter 
(§||), neutrinos (§[|) and the newly-popular vacuum 
energy (§||). In §|7| we consider the hypothesis (ad- 
vanced by one of the authors for ten years) that 
vacuum energy is in fact so dominant that there is 
no need for an exotic component to the dark mat- 
ter. Conclusions are summarized in §||. 

2 Dark Matter and the 

Evolution of the Universe 

The relevant equations here are Einstein's field 
equations of general relativity :Q 

^ 1 . 8vrG _ 

where the terms on the left-hand side describe the 
geometrical structure of spacetime, while those on 
the right-hand side describe its matter and and ra- 
diation content. IZ^ and 1Z are the Ricci tensor 
and curvature scalar, g^ u is the metric tensor, and 
T^ v is the energy-momentum tensor; G, c and A are 

1 Alternative theories of gravity can also be constructed 
in such a way to remove the need for large amounts of dark 
matter |£fj, B5|; we do not consider these here. 
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all constants of nature. About the value of A, in 
particular, we will say more in §|6|. 

The densities of matter (fi M ) and radiation (Q R ) 
are stored in 7^,, while that of the vacuum (Q A ) is 
a function of A.0 Assuming that these three compo- 
nents are distributed homogeneously and isotropi- 
cally on large scales, and that they do not exchange 
energy at a significant rate, Eqs. (|l]) reduce to a 
pair of differential equations in the cosmological 
scale factor R and its time derivatives, including 
the Hubble parameter H = R/R (the expansion 
rate of the Universe) . The equation for H is 

[H(z)/H ] 2 = n M>0 (l + zf + O Rj0 (l + zf + Ov,o 
-(ft TO T,o-l)(l + ^) 2 , (2) 

where z = (R/Rq)^ 1 — 1 is the cosmological red- 
shift. The subscript "0" here (and throughout our 
review) denotes quantities measured at the present 
time; i.e., at redshift z = 0. These are constants, 
and must be carefully distinguished from functions 
of time (or redshift) such as Q M ,f2 R and S7 A . The 
constant f2 TOTj0 = fl M ,o + ^r,o + ^a,o is of particu- 
lar interest, because it separates spatially spherical 
models from hyperbolic ones. If f2 T oT.o > 1, then 
the Universe is closed and finite in extent. Con- 
versely, if f2 TOTj0 < 1, then it is open and infinite 
in extent. And, if TOTj[) = 1 exactly, then we live 
in an infinite Universe which is spatially flat ( "Eu- 
clidean" ) .^| To determine the shape of the homoge- 
neous and isotropic world has been a prime goal of 
cosmologists since expansion was discovered. 

Eq. (||) already tells us a great deal about the 
evolution of the Universe. The first term on the 
right-hand side, O M (l + z) 3 , shows that matter 
acts to increase the expansion rate H(z) as one goes 
to higher z — that is, to slow down the expansion 

2 Densities throughout this review will be written in the 
form of the dimensionless density parameter fl, which is 
just the ratio of physical density (p) to the critical den- 
sity Pcrit(i) = 3H 2 (t)/SivG. This latter quantity depends 
on the Hubble parameter, Eq. ([]), and takes the value 
pcrit.o = 3Ho/8-kG at the present time. If H lies in the 
range 70 — 90 km s _1 Mpc -1 (§^|), then p crit ,o is equivalent 
to between 5.5 and 9.1 protons per cubic meter. 

3 Strictly speaking, it is also possible to obtain flat and hy- 
perbolic solutions which are finite, by suitably "identifying" 
different points and adopting a nontrivial topology |82|; we 
do not pursue this possibility here. 



with time. This is the braking effect of matter's 
gravitational self-attraction. 

The second term, f2 R , (l + -z) 4 , shows that ra- 
diation has the same effect, but with a stronger 
dependence on redshift. This means that, as one 
moves backward in time, photons (and relativistic 
particles) become increasingly important compared 
to pressureless matter. In fact, the dynamics of 
the early Universe (at redshifts above z > 1000) 
must have been completely dominated by them. At 
present, however, the total density O Ri0 of "radia- 
tionlike matter" (including both photons and rel- 
ativistic neutrinos) is several orders of magnitude 
below that of nonrelativistic matter.^] Since we are 
largely concerned in this review with redshifts less 
than ten, we will drop the radiation term in Eq. (|2|) 
from this point onward. 

The third term, f2 Aj0 , is independent of redshift, 
which means that its influence is not diluted with 
time. Vacuum energy will therefore eventually 
come to dominate the dynamics of the Universe in 
any model with A > 0. In the limit t —* oo, in 
fact, the other terms drop out, and we find that 
the density of the vacuum may be expressed as 
^a,o = (Hvo/Hq) 2 , or (since O A , = Ac 2 /3# 2 ) 

Ac 2 = 3Hl , (3) 

where is the limiting value of the Hubble pa- 
rameter as t — > oo (assuming that this latter quan- 
tity exists; i.e., that the Universe does not recol- 
lapse in the future) . This provides a little-discussed 
connection between A (an apparent constant of na- 
ture) and the asymptotic expansion rate (a dynam- 
ical parameter). If A > 0, and if we are living 
at sufficiently late times, then Eq. (|3|) immediately 
predicts that we will measure £l Afi ~ 1. 

The fourth term in Eq. (Q), finally, shows that 
an excess of ri TOT over one (i.e., a positive curva- 
ture) acts to offset the contribution of the first three 
terms to the expansion rate, while a deficit (i.e., a 
negative curvature) enhances them. Open models, 
in other words, expand more quickly at any given 

4 This is inferred, not only from measurements such as 
those of the Cobe satellite, but also from the fact that a 
too-high density of relativistic matter would have slowed ex- 
pansion so much that the Universe could not have lived long 
enough to produce the oldest stars we see. 



3 



redshift z (and therefore last longer) than closed 
ones. This curvature term, however, goes only as 
(1 + z) 2 , which means that its importance drops off 
relative to the matter and radiation terms at early 
times, and becomes negligible compared with that 
of the vacuum term at late ones. 

As recently as the 1980s, many cosmologists were 
persuaded that Eq. (||) could be substantially sim- 
plified, not only by neglecting the second (radia- 
tion) term, but also the third (vacuum) and fourth 
(curvature) terms on the right-hand side. This ap- 
peared reasonable at the time, for four principal 
reasons. First, these terms differ sharply from each 
other (and from the first term) in their dependence 
on redshift z, and the probability that we should 
happen to find ourselves in an era when all four 
terms have similar values would seem a priori very 
remote. By this "Dicke coincidence" argument, it 
was felt that only one term ought to dominate at 
any given time p^] . Second, the vacuum term in 
particular was avoided for historical and theoreti- 
cal reasons (to be discussed in §^). Third, a pe- 
riod of cosmic inflation was widely asserted to have 
driven Q TOT (i) to exactly unity in the early uni- 
verse.0 And finally, this "standard Einstein-de Sit- 
ter" (EdS) model was favored on grounds of sim- 
plicity. These arguments are no longer valid today. 
We are justified in neglecting the radiation term, 
and only the radiation term in Eq. (0), leaving 



H\z) = H 2 [ft M , (l + z) 3 + ft A , 

- (0 TOT , -l)(l + z) 2 ] 



(4) 



This is the modern version of what is usually called 
Friedmann's equation in cosmology. It may be inte- 
grated numerically for the cosmological scale factor 
R(t) as a function of time. 

Several examples are plotted in Fig. 2, includ- 
ing closed models (1 through 5) and one flat (6) 
and open model (7). Model 1, with O M = 0.014 
and Q A0 = 1.08, has been proposed in [79, 30 1 and 



J That this is not necessarily so has been shown by sev- 
eral authors fsil , pq| . A suitable inflationary phase may be 
grafted onto models with f2 T oT,o 7^ 1. The probability of 
finding oneself in such a universe depends, not just on the 
raw amount of inflation, but on factors such as the elapsed 
time and the distribution of initial conditions (such as phase 
transitions ]l8|) which preceded the inflationary epoch. 
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Fig. 2. Evolution of the cosmological scale factor R(t)/Ro as 
a function of time in several models with A > 



will be discussed further in §0. The others all have 
= 0.3, a fi gure widely quoted today for the to- 
tal density of gravitating matter (fQ see, e.g., Q). 
Model 6, with ft Aj0 = 0.7 (known as the ACDM 
model), has been singled out as the newest "stan- 
dard model" of cosmology. Two timescales are plot- 
ted (top and bottom), depending on the present 
value H of Hubble's parameter (§||). 

Along each of the curves, we have marked the 
points where Q M ( Z ) takes on maximum values (A), 
the points of inflection (*), and the points where 
f2 A (z) takes on maximum values and H(z) reaches 
a minimum ( V). One finds that R takes the special 
value R A = 1/y/X at the points (a). The cosmolog- 
ical constant may thus be understood physically (in 
closed models) as the curvature of space at the time 
when the matter density parameter goes through its 
maximum (see [^] for discussion). 

Joining the points of inflection in Fig. 2 are two 
dashed lines marked f2 AjE (for "Einstein limit"). 
One must have O A < S1 AE (a function of the mat- 
ter density O M ) in order for expansion to orig- 
inate in a big bang singularity. Solutions with 
r2 A = f2 AE go over to Einstein's static model as 
t — > 00. When U Afi > Sl AiE , R(t) drops to a nonzero 
minimum and starts to climb again in the past 
direction; these are usually known as Eddington- 
Lemaitre (or "bounce") models. 

The value of f2 AjE can be computed for a given 
model by differentiating Eq. (||) at the points of in- 
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flection (*). This leads to a cubic polynomial which 
must be solved parametrically in general |^0|, |9^| but 
has a little-appreciated direct solution for cases in 
which fi M . ^ 0.5 jL7|]. With the conviction that 
^m,o = 1 now fading in the astronomical commu- 
nity, and most cosmologists calling for values of 
^m,o ~ 0.3, it may be worthwhile to dust off this 
formula again. In modern form it reads 125 ] 



1 — O -I- — O 2/3 

x J 'M,0 * o M "° 



1 — ri M0 + yl — 2f2 M o j 



.N 1/3 



1 ^M,0 -y/l — 2O Mj0 j 



\ 1/3' 



(5) 



One finds that A , E = 1.10 if O M , = 0.014, for 
instance (§0), while fi Mj0 = 0.06 leads to A E = 
1.25. Combined analysis of data on cosmic mi- 
crowave background (CMB) fluctuations from the 
Cobe, Boomerang and Maxima experiments im- 
plies that TOT(0 ^ 1.24 at 2a confidence |3^]. Since 
ri Ai0 is certainly less than f2 T oT,0) we infer that 
f2 A0 < f2 AE in models with O M>0 ^ 0.06, which 
may be regarded as a proof of the existence of the 
big bang in these models [33]. For higher matter 
densities 0.2 ^ O Mi0 ^ 0.5 (as discussed in §||), one 
obtains even larger values of f2 AjE , between 1.5 and 
2.0 (see the Einstein limits in @). 

The differences between the models shown in 
Fig. 2 become apparent when their evolution is plot- 
ted on a phase diagram, with matter density param- 
eter along one axis and vacuum density parameter 
along the other. The key equations are 



n A ( z ) 



^A,o/ [^M,o(l + + ^A,0 

- (^ Mi0 +^ Ai0 -l)(l+z) 2 ] 

{n Mfi /n Afi )n A (z){i + zf . 



(6) 



Fig. 3 depicts the same family of models as Fig. 2, 
with redshift factors [1+z = R /R(t)] labelled at in- 
tervals along the curves. Also marked are contours 
of constant deceleration, defined by q = —RR/R 2 = 
il M /2 — ^ A - This parameter takes values of 0.5 at 
each point of maximum matter density parameter 
(a), zero at the inflection points (*), and —1 at the 
points of minimum expansion rate (v). 

All positive-A models begin in Fig. 3 at the point 
(1,0) and evolve asymptotically toward (0,1) as 




Fig. 3. Evolution in phase space denned by flm(z) and SIa(z) 
of the models whose scale facto rs a re plotted in Fig. 2. Plots 
of this kind go back at least to [112], and have been extended 
into the fin -direction in |33L 



t — * oo. Flat (Euclidean) models follow a straight 
line; any deviation from critical density produces a 
curved path. Those to the right of Model 6 are all 
closed. Models 5 through 2 are increasingly unlikely 
insofar as they violate the above-mentioned obser- 
vational bound r2 TOT0 ^ 1.24 on total density pq| . 
Model 2, in particular, cannot describe the real uni- 
verse. However, the model immediately adjacent to 
it in phase space (Model 1) is perfectly acceptable in 
this regard, since it has T ot,o = 1-094. Very differ- 
ent combinations of fi M ,o and f2 Ai0 , in other words, 
can produce almost identical trajectories in phase 
space. Indeed, from the perspective of Fig. 3, the 
popular Euclidean models appear implausibly fine- 
tuned. One would not expect the Universe to take 
the shortest path through phase space, any more 
than one would expect a star to follow a straight 
line from the main sequence to the red giant branch 
on a Hertzsprung-Russell diagram]^ 

The slow expansion rate and high matter den- 
sity parameter between the points marked (a) and 



6 Of course, the Universe we have assumed (homogeneous 
and isotropic) is simpler than most stars, and might be re- 
stricted to such a path for reasons having to do with some 
higher symmetry of nature. While this may be so, however, 
any such symmetry would lie outside the context of Einstein's 
theory of gravity as it stands. 
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( V) single out this stage of evolution for large-scale 
structure formation. If J7 A is of the same order 
(or less) than fi Mj0 , however, this process must oc- 
cur very quickly. Consider ACDM, represented by 
Model 6 in Fig. 2. Analysis of the Hubble and 
William Herschel Deep Fields suggests that the 
number density of galaxies at redshifts z ~ 4 — 6 
(i.e., at scale factors R/R ~ 0.1 — 0.2) is equal to 
that at z = ||109[| . If so, then these objects were 
in place by z = 4, or (consulting Fig. 2 for Model 6, 
and using either the top or bottom scale for H ) 
within 1.2-1.5 Gyr after the big bang. This poses 
a serious challenge for the model, since primordial 
density fluctuations must not only decouple from 
the Hubble expansion on short timescales, but do 
so at a time when the expansion rate (the slope of 
the curve in Fig. 2) is some six times its present 
value. The problem is even worse in models with 
lower values of O A (e.g., Model 7). 

The standard way to address this has been to sup- 
pose that most of the matter density is in an exotic 
new form (CDM) which is able to decouple from 
the primordial fireball before the baryons, prepar- 
ing potential wells for them to fall into (§|j). This 
approach successfully accelerates structure forma- 
tion on large scales [j94j| , but is perhaps too suc- 
cessful on smaller ones. Galaxy-sized regions are 
formed with excessively peaked central masses (the 
"density cusp problem") [|87| and too many low- 
mass fragments (the "substructure problem") |74| . 
These problems may be resolved within the CDM 
picture by refining the properties of the exotic mat- 
ter; proposals include warm [63], fuzzy |jB5f] , fluid 



SCDM 



ACDM 



pq ] and self- interacting dark matter 1 1 11 

Alternatively, difficulties with the growth of 
large-scale structure are substantially eased in mod- 
els with larger ratios of O A , to U M>0 [39, 106 1. In 
Model 1, for instance, Fig. 2 shows that redshift 
z = 4 corresponds to between 9.3 and 11.9 Gyr af- 
ter the big bang (depending on the value of H ), 
giving the galaxies seven times longer to form. The 
expansion rate at this redshift is only 0.7 times its 
present value. Nor is the low (present) density of 
gravitating matter a problem in this model, because 
f2 M (z) reaches levels as high as four times the crit- 
ical density at redshifts near z « 5 (Fig. 3). It is 
natural to associate this redshift with the onset of 





Fig. 4. Numerical simulations of structure formation. On 
the left is the "old standard model" (SCDM) with Q M ,o ~ 
1, f2 A ,o = and h — 0.5. On the right is the "new standard 
model" (ACDM) with n M ,o = 0.3, fi A ,o = 0.7, and ha = 0.7. 
Panel size is scaled to match the Hubble expansion; time 
runs from top (z = 3) to bottom (z = 0). (Images courtesy 
J. Colberg and the ViRGO Consortium). 



galaxy formation, and it would be of great interest 
to test "slightly-closed" cosmologies of this kind via 
numerical simulations. 

Existing studies like that shown in Fig. 4 by the 
ViRGO Consortium [p^] have so far been restricted 
to flat and open models. Of these, ACDM (the 
new "standard CDM" model) produces a mass dis- 
tribution in much better agreement with the ob- 
served distribution of light than EdS (the old "stan- 
dard CDM" model, still denoted SCDM in many 
works). The improvement is especially pronounced 
at higher redshifts (top two panels). This has been 
taken as evidence for a significant A-term. More 
detailed analysis, however, reveals that the power 
spectrum of the mass distribution does not agree 
with that observed for galaxies in either of these two 



models [39|. The discrepancies have typically been 
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attributed to nonlinear, scale- and morphology- 



dependent bias factors [16|; a complementary ap- 
proach would be to consider an expanded repertoire 
of J7 Mi0 and Q A]0 values. 

3 Baryons 

Beginning from the observed luminosity density of 
the Universe, one can infer the total density of lumi- 
nous baryonic matter (i.e., that in stars) by making 
various reasonable assumptions about the fraction 
of galaxies of different morphological type, their ra- 
tios of disk-type to bulge-type stars, and so on. The 
latest such estimate is [H91 



n hVM = (0.0027 ± O.OOlOfyT 1 



(7) 



where h is the present value of Hubble's parameter 
expressed in units of 100 km s _1 Mpc -1 . 

This latter parameter is unfortunately not yet 
fixed by observation, and we pause to discuss its 
value before proceeding. Using various relative- 
distance methods, all calibrated against Cepheid 
variable stars in the Large Magellanic Cloud 
(LMC), the Hubble Key Project (Hkp) team has 
determined that h = 0.71 ± 0.06 [p8[ . Fundamen- 
tal physics methods (e.g., the Sunyaev-Zeldovich 
effect, gravitational lensing time delays [GLTDs]) 
have higher uncertainties but are roughly consistent 
with this, h = 0.65 ± 0.08 [10C]. The near conver- 
gence of these approaches has been widely hailed, 
with many authors asserting that precision values 
of h are just around the corner, awaiting only a 
final round of experimental refinements. 

On the other hand, a recalibrated LMC Cepheid 
period-luminosity relation based on a much larger 
Cepheid sample (from the Ogle collaboration) 
leads to a considerably higher value of h = 0.85 db 
0.05 [126]. And a new VLBI measurement of the 



transverse velocity of water masers in NGC4258 
gives a purely geometric distance to this galaxy 
[ |6C| l which also implies that the traditional cali- 
bration is off, boosting Cepheid-based estimates by 
12 ± 9% Q. This would raise, e.g., the Hkp value 
to h = 0.80 ± 0.09. There is some additional in- 
dependent support for this recalibration in observa- 
tions of eclipsing binaries and "red clump stars" 



in the LMC 119 [. Fundamental physics approaches 
are also not immune to systematic effects: GLTD- 
based values of h , which are routinely computed 
assuming EdS, rise by about 7% (on average) in 
ACDM, and 9% in open models. 

"Hubble fatigue" may therefore have prompted 
cosmologists to embrace prematurely small levels of 
uncertainty in h . We attempt to allow for this by 
retaining two possible values for h in this review; 
a "low value" of h = 0.7 and a "high value" of 
h = 0.9. This seemingly modest range of choices 
turns out to discriminate quite powerfully between 
the cosmological models considered here. To a large 
extent this is a function of their ages. Fig. 2 reveals, 
for example, that ACDM (represented by Model 6) 
is 13.5 Gyr old if h = 0.7, or 10.5 Gyr if h = 0.9. 
The oldest metal-poor halo stars seen in the Milky 



Way have an age of 15.6±4.6 Gyr [27], setting a firm 
lower limit of 11.0 Gyr on the age of the Universe. 
This is enough to rule out ACDM with the high 
value of h , but not the low one. 

Model 1, on the other hand, faces the opposite 
problem: Fig. 2 shows that it has a total age of 
30.2 Gyr if h = 0.9, or 38.8 Gyr if h = 0.7. 
Both numbers are larger than most cosmologists 
are prepared to accept. However, upper limits on 
the age of the Universe are not as secure as lower 
ones. One must take into account, for instance, 
that the galaxy formation associated above with 
redshifts z ~ 4 occured between 9 and 12 Gyr after 
the big bang in this model, so that galaxies would 
not be older than 24 ± 3 Gyr in any case. 

There are various ways to test such a hypothe- 
sis. One might expect to find a greater spread in 
galaxy ages (and hence colors) at z ~ 4, given their 
longer formation time. Galaxies in Model 1 had 
~6 Gyr to form, or about ~25% of their nominal 
lifetime, according to Fig. 2 (with h = 0.9). The 
corresponding fraction in Model 6 is <10%. This 
may not necessarily translate into a difference be- 
tween observed color spreads, however, since high- 
redshift galaxies are seen principally during (rela- 
tively brief) episodes of star formation. 

Very old galaxies, if they exist, should also be 
present at lower redshifts. They would be inher- 
ently faint and reddened, making them difficult to 
find and distinguish from younger objects which are 
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simply obscured by dust. Nevertheless, several can- 
didates have been noted in the past years, including 
a number of low surface brightness galaxies and 
extremely red objects |25|, [34| whose ages based on 
simple ("single-burst") evolution models appear to 
be as high as ~17 Gyr or more. 

If our own Milky Way is not unusually young, we 
should also expect to find large numbers of dead 
stars in the galactic halo. These would act as mi- 
crolenses, inducing variability in background stars 
and quasars, even if they were too dim to be seen 
directly. Such objects may now have been detected 
in the direction of the LMC (see below). 

Returning now to the density of luminous matter, 
we find with our values for h that Eq. (0) gives 

fi LUM = 0.0034 ± 0.0018 . (8) 

The visible Universe, in other words, constitutes an 
insignificant 0.5% or less of the critical density. 

It may, however, be that many of the baryons 
in the Universe are not visible. How significant 
could these dark baryons be? The theory of cos- 
mic nucleosynthesis provides us with an indepen- 
dent method for determining the density of total 
baryonic matter in the Universe, based on the as- 
sumption that the light elements we see today were 
forged in the furnace of the hot big bang. Results 
using different light elements agree tolerably well, 
which is impressive in itself. The primordial abun- 
dances of 4 He (by mass) and Li (relative to H) im- 
ply a baryon density of f2 BA R = (0.011 ± 0.005)h~ 2 
pQ ], whereas new measurements based exclusively 
on the primordial D/H abundance give a higher 
value: O BA r = (0.019 ± 0.002)/i,7 2 These two 

results, both given at the 2a confidence level, are 
superimposed on a plot of predicted light element 
abundances in Fig. 5. They do not overlap. In our 
view it is premature at present to exclude either 
one. We therefore adopt the same strategy here as 
with Hubble's parameter, retaining a "low" baryon 
density of O.OI/17 2 and a "high" one of 0.02V 2 
throughout our review. Combining this with our 
high and low values of h , we conclude that the 
baryonic density lies in the range 

Q BAR = 0.012 - 0.041 , (9) 

a result in very good agreement with that obtained 
by the entirely independent method of adding up 
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Fig. 5. Measurements of the primordial abundance of light 
elements 4 He, D and 7 Li (scales on left-hand side) determine 
the value of Obah (top). The bottom scale gives the corre- 
sponding baryon-to-photon ratio r\ (adapted from 115]). 



individual contributions from all known repositories 
of baryonic matter via their estimated mass-to-light 
(M/L) ratios |Q. If O TOT:0 is close to unity, as it 
now seems (§|6|), then it follows from Eq. (||) that 
baryons — and everything that would have been 
recognized as "matter" before 1930 — make up less 
than 5% (by mass) of the known Universe. 

Most of these baryons, moreover, have not been 
seen. The baryonic dark matter fraction / BDM (= 
^bdm/^bar) = 1 - ^lum/^bar lies in the range 



/b 



77% 



o , 



(10) 



where we have used Eq. (Q) together with our high 
and low values of h and f2 BAR . Where could these 
dark baryons be? One possibility is that they 
are smoothly distributed in a gaseous intergalactic 
medium, which would have to be strongly ionized 
in order to explain why it has not left a more obvi- 
ous signature in quasar absorption spectra. Recent 
observations using Ovi absorption lines tracer 
of ionization suggest that the contribution of such 
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material to n BAR is at least O.OOS/i^ 1 [ 116 |, com- 
parable to r^LUM- Numerical simulations are able 
to reproduce many observed features of the Lyman 
a (Lya) forest with as many as 80 — 90% of the 
baryons in this form [86]. 

Dark baryonic matter, however, could also be 
bound up in clumps of matter such as substellar 
objects (jupiters, brown dwarfs) or stellar remnants 
(white, red and black dwarfs, neutron stars, black 
holes). The former are not thought to be numer- 
ous enough to be important, given their low mass. 
The latter are limited in the opposite sense; black 
holes cannot be more massive than about 1O 5 M 
since this would lead to dramatic tidal disruptions 
and lensing effects which are not seen [117]. The 
critical mass range for dark baryon clumps is thus 
within a few orders of magnitude of the solar mass. 
Microlensing constraints based on quasar variabil- 
ity do not seriously limit such objects at present, 
setting an upper bound of 0.1 (well above f2 

bar ) tm 

their combined contributions to the density param- 
eter of an EdS universe [107]. 

A likely detection of dark compact objects within 
our own galactic halo has recently been reported 
by the Macho microlensing survey of stars in the 
LMC j|]. The lenses, with masses in the range 
(0.15 — O.9)M , appear to account for between 
8% and 50% of the high rotation velocities seen 
in the outer parts of the Milky Way — depend- 
ing on the halo model chosen, and extrapolat- 
ing (at 2a confidence) from the ~15 events actu- 
ally seen.[] The identity of these objects has been 
hotly debated, with some authors linking them to 
faint, fast-moving objects apparently detected in 
the Hubble Deep Field [67]. It is unlikely that 



they could be traditional white dwarfs, since these 
are formed from massive progenitors whose metal- 



rich eject a we do not see [41]. Degenerate "beige 
dwarfs," which might be able to form above the 
hydrogen-burning mass limit of 0.08 Mq without 
fusing ^] are one possibility. Another would be 
a population of ancient, low-mass (< O.6M ) stars 
which have simply cooled into invisibility^] 



4 Exotic Dark Matter 

Three main reasons have been proposed for go- 
ing beyond dark baryons and introducing a sec- 
ond species of invisible matter, the exotic cold dark 
matter (CDM), into the Universe: (1) a range of 
observational arguments imply that the density pa- 
rameter of total gravitating matter (f2 M ,o = 0, BAR + 
+ is higher than that provided by baryons 
and neutrinos alone; (2) our current understand- 
ing of large-scale structure formation requires the 
process to be helped along by quantities of cold 
(i.e., nonrelativistic) gravitating matter in the early 
universe, creating the potential wells for infalling 
baryons; and (3) theoretical physics supplies sev- 
eral plausible (albeit still hypothetical) candidate 
CDM particles with the right properties. 

Since our ideas on structure formation may yet 
change, and the candidate particles may not mate- 
rialize, the case for exotic CDM turns on the ob- 
servational arguments. At present these agree to 
within no better than a factor of five, pointing to 
values of O M between about 0.1 and 0.5. (Not long 
ago, in the 1980s, there were calls for f2 M ,o = L hut 
these tended to come from theorists wishing to re- 
tain the EdS model, and are no longer tenable ob- 
servationally [125]-) The lower limit is crucial: if 

^M,0 > ^BAR. + ^I>) then f^CDM > 0. 

The arguments can be broken into two classes: 
those which are purely empirical, and those which 
assume in addition the validity of the gravitational 
instability (GI) picture of structure formation. Let 
us begin with the empirical arguments. One has 
been encountered already: the galactic rotation 
curve. If the Macho results are taken at face value, 
and if the Milky Way is typical, then it is probable 
that dark compact objects make up less than 50% 
of the halo mass in spiral galaxies. The remaining 
halo dark matter does not appear to consist of bary- 
onic matter in known forms such as dust, rocks, hot 
or cold gas, or hydrogen snowballs |j59| . 

The total amount of dark matter in spiral galax- 
ies, however, is rather limited. The easiest way 



7 An alternative interpretation of the data, that most or 
all of these lenses are actually faint stars inside the LMC 
itself, now appears to be strongly disfavored |89| . 

8 Existing limits on the density of halo objects in this mass 



range |>2j necessarily extrapolate in a number of ways from 
known stellar populations, and are based on theoretical stel- 
lar evolution models J|, which do not currently extend 
beyond 20 Gyr. 
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to see this is to divide the total dynamical mass 
of the Milky Way (including its unseen halo mat- 
ter) by its luminosity.[] The resulting mass-to- 
light ratio, M/L = (21 ± 7)M /L , is a mere 
~ Q.02/1^" 1 times that of a critical-density universe, 



(M/L) ciit = (1136 ± 138)/i o M /L ]2C 

Most of the mass of the Universe, in other words, 
is spread over scales larger than galaxies, and it is 
here that the arguments for exotic CDM are most 
compelling. The M / L-ratio method is in fact the 
most straightforward: one measures M/L for a cho- 
sen region, corrects for the corresponding value in 
the "field," and divides by (M/L) crit to obtain f2 M ,o- 
Much, however, depends on the region. A widely 
respected application of this approach, that of the 
Cnoc team, uses rich clusters of galaxies. These 
systems sample large volumes of the early Universe, 
have dynamical masses which can be measured by 
three independent methods (the virial theorem, x- 
ray gas temperatures and gravitational lensing), 
and are subject to fairly well-understood evolution- 
ary effects. They are found to have field M/L-ratios 
of (213 ± 59)/i o M /L , giving O M , = 0.19 ± 0.06 
(1<7 confidence) when Q Afi = [2C]. This result 
scales as (1 — 0.4f2 Aj0 ) @| so that, e.g., S7 M ,o drops 
to 0.11 ± 0.04 when Q Afi = 1. 

The weak link in this chain of inference is that 
rich clusters may not be characteristic of the Uni- 
verse as a whole. Only about 10% of galaxies are 
found in such clusters. If individual galaxies (like 
the Milky Way, with M/L « 21) are substituted for 
clusters, then the inferred value of £l M ,o drops by a 
factor of ten, approaching BAR (§[|) and removing 
the need for exotic CDM. A recent comprehensive 
effort to address the impact of scale on M/L ar- 
guments concludes that Q Mfi = 0.16 ± 0.05 when 
regions of all scales (from individual galaxies to su- 
perclusters) are considered 

A second line of argument uses the cluster baryon 
fraction (M BAR /M TOT ) of baryonic to total gravitat- 
ing mass in clusters. Baryonic matter is defined as 



9 M MW has been calculated at (4.9 ± 1.1) x 10 n Mo inside 



50 kpc using using the motions of galactic satellites |75j ; while 
^mw is given by (2.3 ± 0.6) x 1O 1O L in the B-band. 

10 This result is based on comparison with simulations that 
assume Q M ,o = 0.37 plus spatial flatness. It is hence subject 
to systematic uncertainty (in the downward direction) of as 
much as -20% lid], possibly more if Qa,o > 0.63. 



the sum of visible galaxies and hot cluster gas (the 
mass of which can be inferred from its x-ray tem- 
perature). Total cluster mass is measured by one 
or all of the three methods listed above (virial, ar- 
ray, or lensing). At sufficiently large radii the clus- 
ter may be taken as representative of the Universe 
as a whole, so that Q M ,o = BAR /(M BAR /M TOT ), 
where n BAR is fixed by nucleosynthesis (§|3|). Ap- 
plied to various clusters, this procedure leads to 
^m,o = 0.3 ±0.1 || — a result which is almost cer- 
tainly an upper limit, partly because baryon enrich- 
ment is more likely to take place inside the cluster 
than out, and partly because dark baryonic mat- 
ter (e.g., Machos) is not taken into account; this 
would raise M BAR and lower f2 M ,o- 

A final, recent entry into the list of purely em- 
pirical methods uses the separation of radio galaxy 
lobes as standard rulers, a variation on the clas- 
sical angular-size distance test in cosmology. The 
widths, propagation velocities and magnetic field 
strengths of these lobes have been measured for 
14 radio galaxies with the aid of long-baseline ra- 
dio interferometry, leading to the (la) constraint 
fi M ,o < 0.10 for n Afi = 0, or ft M>0 = 0.10+^ for 
flat models (fi A ,o = 1 - fi M ,o) PRl- 

We consider next the Gl-based measurements of 
r2 M , which are "circular" in the sense that they as- 
sume that large-scale structure formed via gravita- 
tional instability from a Gaussian spectrum of pri- 
mordial density fluctuations — a process which (as 
we currently understand it) could not have taken 
place as it did unless £l M .o is considerably larger 
than f2 B AR- According to GI theory, formation of 
large-scale structure is more or less complete by 
z ai tt^o - 2 Hi- Therefore one can constrain f^M.o 
by looking for evidence of number density evolution 
of large-scale structures such as galaxy clusters. In 
a low-density universe, this would be relatively con- 
stant out to at least z ~ 1, whereas in a high-density 
universe one would expect the abundance of clus- 
ters to drop rapidly with z because they are still in 
the process of forming. In fact, massive clusters are 
seen at redshifts as high as z = 0.83, leading to the 
(la) limits f2 M>0 = 0-17io*09 f° r ^a,o = models, 
and ft M ,o = 0.2218:5? for flat ones |]. 

Evolution of the mass power spectrum P(k) con- 
strains £7 M o in a similar way. Here one uses the 
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fact that structures of a given mass form by the 
collapse of large-scale regions in a low-density uni- 
verse, or smaller-scale regions in a high-density one. 
Comparing P(k) for the present-day distribution of 
galaxy clusters to that for the distribution of matter 
at some earlier epoch therefore yields an estimate 
of f2 M ,o- Using the mass power spectrum of Lya 
absorbers at z ~ 2.5, for instance, one finds that 
Q M ,o = 0.46±g;ig {la) for ft A:0 = models pj| . 
This result goes as approximately (1 — 0.417^,0)) so 
that the central value of O M drops to 0.34 in a flat 
model, and 0.28 if Q Afi = 1. 

A final group of measurements comes from galaxy 
peculiar velocities. These are produced by the grav- 
itational potential of locally over- (or under-) dense 
regions relative to ft M ,o, but also depend on fi Mj0 
itself. The power spectra of the velocity and den- 
sity distributions can be related within the context 
of GI theory. A typical bound derived from sev- 
eral such studies is Cl M .o > 0.3 |128|. Dependence 



on f2 A is modest since these tests probe relatively 
small volumes, but lower limits derived in this way 
can depend significantly on h as well as the spectral 



index n of the density distribution. In [125], where 
the latter is normalized to CMB fluctuations, re- 
sults take the form tt^hl^n 2 m 0.33 ± 0.07 (2a). 
The preferred value of Q M ,o therefore drops from 
0.53 ± 0.11 (if h = 0.7) to 0.38 ± 0.08 (if h = 0.9), 
where we have assumed n = 1. 

To summarize, one may say that purely empir- 
ical arguments lean toward values of £l M ,o ~ 0.3 
or lower, whereas GI (gravitational instability) 
theory-based results tend to come in at ~ 0.2 and 
higher. If there is flexibility in the lower limits on 
O M , it lies in the empirical methods, especially if 
contributions from dark baryons are near their up- 
per limit (§||). It is unlikely, however, that the limits 
based on GI theory can be stretched far enough to 
remove the need for exotic CDM. We therefore con- 
clude that this component of the dark matter has a 
density parameter in the range 



0.1 

- 



- 0.5 
0.4 



(GI theory) 
(otherwise) 



(11) 



If our current understanding of structure formation 
via gravitational instability is correct, then exotic 
CDM must exist. Conversely, if exotic CDM does 



not exist, then our understanding of structure for- 
mation is incomplete. 

The debate, of course, becomes moot if exotic 
CDM (with f2 C DM ~ 0.3) is actually discovered in 
the laboratory. Theorists have proposed a colorful 
list of particle candidates, with varying degrees of 
testability. Two have emerged as most plausible: 
the axion and the weakly interacting massive parti- 
cle (WIMP)] 11 ] Either one of these could in princi- 
ple make up the CDM because each, if it exists, is 
(1) cold (i.e., nonrelativistic in the early Universe), 
and (2) expected on theoretical grounds to have a 
collective density within a few orders of magnitude 
of the critical density (see Appendix). Ambitious 
experimental detection efforts around the world are 
now directed at both particles. While they have not 
turned up anything so far, most of the theoretical 
parameter space remains unexplored. 

5 Neutrinos 

Since neutrinos indisputably exist, and in great 
numbers (their number density n u is 3/11 times 
that of the CMB photons, or 112 cm -3 per species), 
they have been leading particle dark matter candi- 
dates for longer than either the axion or the WIMP. 
They gained prominence in 1980 when teams in 
the U.S.A. and Soviet Union both claimed to have 
evidence of nonzero neutrino rest masses. While 
these claims did not stand up, a new round of more 
sophisticated experiments is once again suggesting 
that m u (and hence £l u ) > 0. 

Dividing n v m v by the critical density of the Uni- 
verse, one obtains immediately 



(12) 



where the sum is over the three neutrino species -P| 
Current laboratory upper bounds on neutrino rest 

11 Other contenders include WIMPzillas, primordial black 
holes, magnetic monopoles, solitons (in various dimensions), 
cosmic string loops, shadow matter and mirror matter. We 
refer the interested reader to the latest confe rence reports 
|t2^ [r^ l or semipopular works such as J7f], 
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The calculations in this section are strictly valid only for 
m v c 2 < MeV. More massive neutrinos with m„c 2 ~ GeV 
were once considered as CDM candidates but are no longer 
viable since LEP experiments rule out additional neutrino 
species with m„c 2 < 46 GeV (i.e., half the Z -mass). 
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masses are 3 eV (m^c 2 ), 0.19 MeV (m^c 2 ) and 
18 MeV (m UT c 2 ), so it would appear feasible in prin- 
ciple for these particles to close the Universe. In 
fact m„ and m Vr are limited far more stringently 



by Eq. (12) than by laboratory bounds. 

Perhaps the best-known theory along these lines 



in recent years is that of Sciama jlOql , who postu- 
lated a population of 29 eV r-neutrinos, decaying 
via v T — > Vn + 7 into much lighter /i-neutrinos and 
15 eV photons on timescales of t v < 3 x 10 23 s. 
Decay photons with these properties would solve 
a number of astrophysical puzzles, such as the 
high degree of ionization in the interstellar medium 
and the large intergalactic flux of hydrogen-ionizing 
photons. The proposed neutrinos would moreover 
provide exactly the critical density if h Q = 0.56. 
This model, however, has now been ruled out by 
the absence of a strong 15 eV emission line in the 



extragalactic background light |92|| . 

More generally, the strongest upper limits on Q, v 
come from our current understanding of structure 
formation. Neutrinos are hot (i.e., relativistic at 
the time of decoupling from the primordial fireball) 
and therefore able to stream freely out of density 
perturbations in the early Universe, erasing them 
before they have a chance to grow. Good agree- 
ment with observations of large-scale structure can 
be achieved in models with Q u as high as 0.2, but 
only if n BAR + f2 CDM = 0.8 and h = 0.5 pjf. A 
more realistic upper limit follows from a statistical 
exploration of the entire parameter space and leads 
to the conclusion that m u c 2 < (9.2 eV)f2 CDM (2a) 
over ft CDM ^ 0.6 for flat models ||. Eq. © 
then implies that £l u < 0.12S1 CDM (if h = 0.9) or 
0.20S1 CDM (if h = 0.7) — a neutrino density below 
that attributed to exotic CDM (§§), but still well 
above that of the baryons (§^).0 

Unexpected new lower limits on £l u have come 
from atmospheric (Super-Kamiokande |44|]), solar 
(Sage 0], Homestake pi, Gallex Q), and 
accelerator-based (Lsnd ||) neutrino experiments. 



13 These (possibly quite significant) values of Q„ do not 
invalidate our earlier assumption (§^|) that the present den- 
sity of radiationlike matter is negligible. While neutrinos 
are relativistic at decoupling, they lose energy and become 
nonrelativistic on timescales t NR ~ 190, 000(m„c 2 /eV)~ 2 yr 
fr^ | — well before the present epoch for neutrinos which are 
massive enough to be of interest. 



In each case it appears that two neutrino species 
are oscillating into each other, a process which can 
only take place if m u > 0. The strongest evi- 
dence comes from Super-Kamiokande, which has 
reported oscillations between r- and ^-neutrinosp^ 
with 5 x 10" 4 eV 2 < Am^c 2 < 6 x 10~ 3 eV 2 (2a), 
where Am^ = |m 2 r — m 2 J [44]. If neutrino masses 
are hierarchical, like the masses of other fermions, 
then ui Ut ^> m„ and m UT c 2 > 0.02 eV. In this 
case it follows from Eq. (jT|) that Q u > 0.0003 (if 
h = 0.9) or 0.0005 (if h = 0.7). If, instead, neu- 
trino masses are nearly degenerate, then Q u cannot 
be determined from this result, but will in any case 
still lie below the upper bound imposed by struc- 
ture formation above. We conclude that the possi- 
ble range of values for this parameter is 



Q u = 0.0003 - 0.2O- 



(13) 



The neutrino contribution to n TOT0 is anywhere 
from an order of magnitude below that of the visible 
stars and galaxies (§^) up to as much as half that 
attributed to exotic CDM (§||). If Q C dm is small (or 
zero), however, then £l u must lie near the lower, 
rather than the upper end of this range; since a 
large density of neutrinos will interfere with struc- 
ture formation (as noted above) unless something 
like exotic CDM is present to help hold primordial 
density perturbations together. 



6 Vacuum Energy 

There are at least four good arguments for the cos- 
mological constant. The first is mathematical: A 
plays a role in Eqs. ([A) similar to that of the addi- 



tive constant in an indefinite integral [103]. The 
second is dimensional: A specifies the curvature 
radius R A = of a (closed) Universe at the 

moment when the matter density parameter £l M (z) 
passes through its maximum (§|2|), thereby provid- 
ing a fundamental length scale for cosmology (see 
p8| for discussion). The third is dynamical: A de- 
termines the asymptotic expansion rate of the Uni- 
verse according to Eq. (H), Ac 2 = SH^. And the 



14 A second possibility, that of oscillations between and 
a sterile fourth neutrino species (v 3 ), has now been excluded 



at the 3cr confidence level 
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Table 1. Theoretical Estimates of £1 A 



Theory 


Estimated value of p A 


"A,0 


QCD 
EW 
GUTs 


(0.3 GeV) 4 /r 3 c- 5 
(200 GeV) 4 ft- 3 c- 5 
(10 19 GeV) 4 ft- 3 c" 5 


10 16 g cm" 3 
10 26 g cm" 3 
10 93 g cm" 3 


10 44 /i " 2 

W 55 ho 2 

10 122 h~ 2 



fourth is material: A is related to the energy den- 
sity of the vacuum via p A c 2 = Ac 4 /8ttG. 

With all these reasons to take this term seriously, 
why have most cosmologists since Einstein set A = 
0? Computational convenience is one explanation. 
Another is the smallness of most effects associated 
with the A-term. Einstein himself set A = in 
1931 "aus Griinden der logischen Okonomie" — for 
reasons of logical economy — because he saw no 
hope of measuring this quantity experimentally at 
the time. He is often quoted as adding that its in- 
troduction in 1915 was the "biggest blunder" of his 
life ("die grofite Eselei in meinem Leben"). This 
statement, which does not appear anywhere in Ein- 
stein's writings but was rather attributed to him 
by Gamow j|7|], is sometimes interpreted as a re- 
jection of the very idea of a cosmological constant. 
It more likely represents Einstein's rueful recogni- 
tion that, by invoking the A-term solely to obtain 
a static solution of the field equations, he had nar- 
rowly missed what would surely have been one of 
the greatest triumphs of his life: the prediction of 
cosmic expansion. 

The relation between A and the energy density of 
the vacuum has led to a new quandary in more re- 
cent times: the fact that p A as estimated in the 
context of quantum field theories such as quan- 
tum chromo dynamics (QCD), electroweak (EW) 
and grand unified theories (GUTs) implies impos- 
sibly large values of £l A ,o (Table 1). These theo- 
ries have been successful in the microscopic realm. 
Here, however, they are in gross disagreement with 
the observed facts of the macroscopic world, which 
tell us that £l A ,o cannot be much larger than or- 
der unity. This "cosmological constant problem," 
is undoubtedly another reason why many cosmolo- 
gists have preferred to set A = 0, rather than deal 
with a parameter whose microphysical origins are 



while nowhere near the size suggested by Table 1, 
is nevertheless greater than zero. The cosmological 
constant problem has therefore become more baf- 
fling, in that any quantum field-theoretic account of 
this parameter must apparently contain a cancella- 
tion mechanism which is not only good to some 44 
(or 122) decimal places, but which begins to fail 
at precisely the 45th (or 123rd). One possibility 
is to treat A dynamical quantity rather than 
a constant of nature, in which case the observed 
smallness of £l Afi might be attributed to the age of 
the Universe |HJ. It is equivalent to introduce a 
fifth element (known as quintessence) into cosmol- 
ogy [^] . In general, however, this means extending 
Einstein's equations ((!]) to incorporate new (and 
so far unobserved) phenomena such as scalar fields, 
which in turn introduce new terms (intermediate to 
the Q M ,<r and f2 Aj0 -terms) into Eq. (||). While these 
ideas hold theoretical promise, it is likely that they 
too must involve fine-tuning if they are to repro- 
duce exactly the values of Q Afi observed. As an 
alternative explanation, it has been proposed that 
a universe in which fi Aj0 was too large (or small) 
might be incapable of giving rise to intelligent ob- 
servers, so that the fact of our own existence already 
"requires" Q A0 



1 pi . 



still unclear (see [22] for review). 

This, however, is no longer an appropriate re- 
sponse because observations now indicate that fl A , , 



Let us pass now to the observational arguments 
for n A . Some have been mentioned already. It was 
noted in §|| (Fig. 4) that numerical simulations of 
large-scale structure formation match models with 
^a,o = 0.7 better than those with £l A}0 = 0. Ad- 
ditional simulations would allow us to explore the 
parameter space with higher resolution. 

Some of the arguments for exotic CDM in §[| also 
show a (more modest) dependence on O A , - The 
trend in most cases is toward higher values of O A 
in conjunction with lower values of f2 M ,o- in the 
arguments from cluster M/L-ratios and mass power 
spectra, for example, we have seen that raising O A 
from zero to one corresponds to a drop of ~ 40% in 
the preferred value of £l M ,o- 

Tentative lower limits on O A have come from 
galaxy number counts. The comoving volume is en- 
hanced at large redshifts {z > 2) for high-A models, 
leading to greater (projected) number densities at 
faint magnitudes. In practice, it has proven diffi- 
cult to disentangle this effect from galaxy evolution. 
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Fig. 6. Two-dimensional slices in time of the flat EdS model 
(Qm,o = 1, ^a.o = 0) with ho = 0.5 (left), and a closed, low- 
fi M ,o, high-n A , model @ g with h = 0.9 (right). There 
is nearly ten times as much linear distance between redshifts 
z = 3 and 4 in the closed model as there is in the flat one. 
(Both figures to same scale; 1 Gly= 10 9 light-years.) 



Early claims of a best fit at f2 A0 ~ 0.9 [45| have 
been disputed on the basis that the steep increase 
seen in numbers of blue galaxies is not matched in 
the K-band p8j . Attempts to account for evolution 
in a comprehensive way have recently produced the 
2cr lower limit f2 A-0 > 0.53 [114 ], with a good fit (for 
flat models) at Q A 



115 



&Ai0 ~ 0.8 
of constraints across the Cl M 



A preliminary plot 
U A plane |105| ] is 
consistent with these results. 

Measurements of O A>0 from gravitational lens 
statistics are based on a similar premise: the in- 
crease in path length to a given redshift in high-A 
models should mean that more lensed sources are 
seen. Somewhat surprisingly, comparison of the ob- 
served frequency of lensed quasars to that expected 
produces results in poor agreement with those in- 
ferred from galaxy counts; leading in fact to the 
strongest current upper limit on vacuum density: 
O A , < 0.66 (2a) for flat models |7fJ. Dust could 
obscure the distant lenses and allow for a larger 
value of f2 A , ||S5|1 . This objection can however be 
met by moving to radio lenses, which give an only 
slightly weaker 2a bound: Q Afi < 0.73 (for flat mod- 
els) or f2 A>0 < 0.4 + 1.5O M>0 (for nonflat ones) |38| . 
Other sources of potential systematic error remain, 
including source redshift distributions, survey in- 
completeness, lens modelling and evolution in the 
mass profiles of the lensing objects. A recent 2a 
limit from radio lenses is S7 Aj0 < 0.95 pf 
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Lensing provides a second constraint on fi Aj0 in 
closed models, based on the redshift of the antipodes 
(the set of points at radial coordinate distance ir). 



It has been shown [51| that the lensing cross sec- 
tion blows up for sources at this distance, implying 
that the antipodes cannot be nearer than the far- 
thest normally lensed object, currently a galaxy at 
z = 4.92 p^] . It is straightforward to compute the 
antipodal distance in terms of f2 M ,o and ^A.oj one 
finds in this way that f2 A . < 1.57 if f2 M ,o = 0.3, an 
upper limit that drops to < 1.30 if f2 M ,o = 0-1 .P| 
The upper limit is f2 A < 1.10 if O Mi0 = 0.014, es- 
sentially the same as the Einstein limit discussed 
in §0. Fig. 6 shows a two-dimensional slice of this 
model, along with the EdS model drawn to the same 
scale. The antipodes (in the closed case) are at 
z ~ 12, well beyond the highest-redshift objects 
seen to date. Fig. 6 also reveals nearly ten times as 
much linear distance between redshifts z = 3 and 
z = 4 in the spherical model as in the flat one; this 
is why one expects to see more lensed sources (and 
faint galaxies) in a high-A universe. 



15 A recent claim []37[ that this technique "rules out values 
of O a ,o between f .0 and 1.5 for ^ £2m,o ^ 0.3" is therefore 
incorrect. The appropriate upper limit (for ma tter densities 
in this range) is f2 A ,o < 1.5, as reported in Jl23J. 
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The best constraints on fi A come from Type la 
supernovae (SNe la), used as standard candles in 
the classical magnitude-redshift relation, can mea- 
sure A ,o more robustly than either galaxy counts 
or lensing statistics because evolution is less of a 
concern for these objects. The Supernova Cos- 
mology Project |Q an d High-z Supernova Search 
teams |102| 1 have lately caused a great stir with their 
reports that a systematic dimming of SNe la at 
z ~ 0.5 by about 0.25 mag relative to that expected 
in an EdS model is best explained by 

^a,o ~ | ^m,o + l ± \ (2a) . (14) 

This leads at once to a conclusion that most cosmol- 
ogists^j have long been reluctant to consider: that 
we live in a vacuum- dominated universe. Several 
words of caution are in order, however. Observa- 
tions must reach z ~ 2 before one can truly be 
certain of discriminating between models like those 
shown in Fig. 2. Intergalactic dust could mimic 
the effects of a cosmological constant, if it were 
"sifted" during the process of ejection from galaxies 
[||. The neglect of evolution may be more serious 



than claimed [31|. And the physics of SNe la ex- 



plosions needs to be better understood |61]. 

The best constraints on S7 A .o come from CMB 
fluctuations] and in particular, from the new detec- 
tion of the first peak in their angular power spec- 



trum by the Boomerang |2g] and Maxima |57] 
experiments. The location of this peak is a direct 
measure of the largest size of fluctuations in the pri- 
mordial plasma at the moment of last scattering, as 
seen through the "lens" of a curved Universe. The 
two detections, combined with earlier data from the 



Cobe satellite jl4|, imply that j68| 



n A 



l.n 



+0.13 
-0.12 



(2a) 



(15) 



This result is more reliable than all the others dis- 
cussed so far because it bypasses "local" phenom- 
ena such as supernovae, galaxies, and even lensed 
quasars; taking us directly back to the radiation- 
dominated era when physics was very simple. 

Let us therefore use Eq. (15) to calculate the en- 
ergy density of the vacuum. Summing the baryon, 



exotic CDM, and neutrino densities — Eqs. @, 
( |TT|) and (|l^) respectively — gives the total matter 
density f2 M ,o- Substituting this into Eq. (]15|), we 
find the following ranges of values: 



9. 



A,0 



0.3 - 1.1 (GI theory) 
0.4 - 1.2 (otherwise) 



(16) 



Vacuum energy, the new and invisible "fire" of mod- 
ern cosmology, is thus indeed dominant, making 
up the bulk of a dark Universe in which light and 
baryons — the constituents of our familiar world — 
appear almost incidental. 

7 How much exotic dark matter? 

The basis for ACDM as the new favorite among cos- 
mological models lies in the approximate orthogo- 
nality of the CMB and supernova bounds, Eqs. ( |l~4|) 
and (15). Indeed, if we take both of these re- 



sults at face value, we can substitute one into the 
other and solve to find O 



A.O 



0.78 ± 0.23 and 
^m,o = 0.33 ± 0.22. The fact that this latter num- 
ber is very near the center of the range of allowed 
values for f2 M>0 in Eq. (|TT|) has then been taken as 
a further sign of the basic correctness of both the 
ACDM model in particular and the GI theory of 
structure formation in general. 

While this is a self-consistent account, and one 
that agrees with most observations, it suffers from 
one flaw: it is inherently improbable. The densities 
of baryonic matter (and exotic CDM, should it ex- 
ist) evolve at a very different rate from neutrinos; 
and both of these components evolve at very dif- 
ferent rates from vacuum energy. So one has three 
kinds of matter which should not have anything like 
the same density parameters at any given time - 
and yet two of them (at least) do. In the ACDM 
picture, in particular, it seems that we happen to 
live at a time when f2 A ,o and ri M ,o are separated by a 
mere factor of two. To illustrate the unlikelihood of 



3 With some exceptions 



79 



§ 



this "preposterous universe," Carroll [22[] has plot- 
ted the evolution of M an d 0, A for 35 powers of 
ten in scale factor in the past and future directions, 
showing that the probability of finding oneself at 
the moment when they should be within even one 
order of magnitude of each other is exceedingly re- 
mote. 
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It may be misleading to characterize this problem 
in logarithmic terms. Several Gyr were necessary to 
form the first galaxies and stars; and they will all 
be gone after a hundred — so it is natural that 
we should find ourselves within this span of cosmic 
history, at least. In Fig. 7, we have plotted the 
evolution of fJ A (i) and f] M (f) on a linear scale for 
the first ~ 80 (100) Gyr after the big bang in the 
ACDM model with h = 0.9 (0.7). This plot con- 
firms that ACDM is improbable, in the sense that 
the values of Q M and f2 A observed at the present 
epoch (0.3 and 0.7 respectively) are atypical. 

Fig. 7 also shows the evolution of Q A and fi M 
in the vacuum-dominated, A+baryon model dis- 
cussed in several places above (Model 1 in Figs. 2 
and 3), with f2 Aj0 = 1.08 and no exotic CDM 
(f2 M ,o = 0.014 ~ n BAR ). This universe, which was 
originally proposed in (7^, SO | and which we shall 
term here the "Abar" model, is a good deal less pre- 
posterous than ACDM in the sense that a factor of 
~ 80 (rather than two) separates the presently ob- 
served values of vacuum and energy density. Indeed 
these parameters are much closer to their "cosmo- 
logical average" values of one and zero respectively. 
While this in itself does not constitute a case for the 
model, it prompts us to wonder whether A might 
not be more important than most cosmologists have 
been willing to consider. Could vacuum energy be 
not just the dominant, but the only significant com- 
ponent of the dark matter? 

Such an idea would have been unthinkable only a 
few years ago, when it was still routine to set A = 
and cosmologists had two main choices: the "one 
true faith" (O M = 1), or the "reformed" (with each 
believer being free to choose his or her own value 
near O M ~ 0.3). All this has been irrevocably 
altered by the CMB experiments (§§|). If there is a 
single guiding principle in choosing models now, it 
is no longer J7 M « 0.3, and certainly not Q A = 0; 
it is ft M .o + ^a,o ~ 1 from the power spectrum of the 
CMB. With this in mind, we will devote this final 
section of our review to exploring the feasibility of 
the Abar model, which has f2 M ,o + ^a,o = 1-094, in 



TIME t/10 9 yrs (H = 70 km s" 1 Mpc" 1 ) 



excellent agreement with Eq. (15). 

To begin with, any model of this kind, with a 
small £l M ,a and a high S7 Ai0 , must face three main 
objections: (1) the lower limits on Q M in §E[ (2) the 
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Fig. 7. Evolution of f2 M and Qa in the Abar and ACDM 
models (Models 1 and 6 in Figs. 2 and 3). Time is set to 
zero at the present epoch; £bb , the time of the big bang, 
is calculated using h = 0.9 for Abar (bottom scale) and 
h = 0.7 for ACDM (top scale). Compare Fig. 3. 



upper limits on f2 A in §||; and (3) the age problem 
in §||. Let us briefly review these points before mov- 
ing on to what are, to us, the strongest arguments 
in favor of the Abar idea. 

We have argued in §[| that the lower limits on 
r2 M are of two kinds: those which are framed in 
the context of gravitational instability theory (thus 
tacitly assuming the existence of large amounts of 
CDM) and those which are not. Of these, the for- 
mer are certainly incompatible with M = 0.014. 
The latter, however, are more flexible, especially 
when their f2 Aj0 - and /^-dependence is taken into 
account (as it often is not). The distinction is im- 
portant, because CDM may not be needed in the 
Abar model: structures have far longer to form, 
and they do so at a time when expansion is slower 
and densities are higher. 

The observational constraints on f2 Aj0 in §|6| are 
not much more secure than those on Q M ,o- Super- 
novae favor a value of O Aj0 closer to ~ 0.8, but 
have large (and possibly underestimated) uncer- 
tainty factors. The antipodal redshift argument 
restricts us to f2 Aj0 < 1.10, which is just above 
the value we have adopted — suggesting that this 
could become a useful test of the theory when ob- 
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servations are eventually pushed to z « 12 (Fig. 6). 
Gravitational lensing statistics remain the strongest 
argument against values of f2 A as large as that 
considered here. It would be a more compelling 
one, however, if we understood why galaxy counts 
(which rely on essentially the same reasoning) lead 
to such different conclusions. 

Finally, as we have argued in §|2| and §||, the age of 
the Abar universe should be seen as an asset rather 
than a liability, particularly in the area of structure 
formation. One also expects to find a faint popula- 
tion of very old, dead stars; these may be the halo 
objects whose existence is implied by the Macho 
survey. The most direct way to rule out the Abar 
model based solely on age considerations would be 
to prove that h is less than or equal to what we 
have referred to as the "low" value (0.7) above. Al- 
ternatively, if evidence continues to mount for the 
"high" value (0.9), then the "age problem" will be- 
gin to put pressure on the ACDM, rather than the 
the Abar model. 

Let us turn now to the observational case for the 
Abar idea, which rests on two main lines of evi- 
dence: Lya absorption spectra, and the lack of a 
second peak in the power spectrum of CMB fluctu- 
ations. We consider these in turn. 

The forest of Lya absorption lines was first used 



to measure £l M ,o and f2 A]0 in [79, 8C]; the most re- 
cent review of this method appears in |120| 1 . The 
idea is extremely simple: one supposes that Lya 
absorbers, like galaxies, are distributed with a cell- 
like structure, and that absorption lines are pro- 
duced when the line of sight to a distant quasar 
cuts through the cell walls (Fig. 8). The crucial as- 
sumption is then made that the cells expand with 
the Hubble flow, and that evolution within them is 
secondary and can be neglected. This is not un- 
reasonable, given that the expansion velocity of a 
typical cell would be of order ~ 3000 km/s, whereas 
peculiar motions inside the cell walls might be no 
more than about ~ 300 km/s. This picture is also 
in accord with the latest observational work on Lya 
absorbers which, although still tentative, suggests 
that they are indeed distributed in structures which 
expand with decreasing redshift pC|], and have co- 
moving size ~ 26/1" 1 Mpc at z = 2.6 fll^ f 7 ! 

17 In fact, neither assumption may be strictly necessary; a 
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Fig. 8. Schematic illustration of a network of cell-like struc- 
tures, with Lya absorption lines seen by an observer in the 
spectrum of a distant quasar due to the fact that the line of 
sight passes through the cell walls. 

One then simply counts the absorption lines and 
measures the mean spacing AA between them. This 
gives the redshift spacing Az(z) of the cells as a 
function of redshift, which may in turn be related 
to their comoving coordinate size A% by 



[Az(z)f 



H 2 (z) 



(17) 



where \ ls the radial coordinate distance and H 2 (z) 
is given as a cubic in (1 + z) by Eq. (||). Since 
Az(z) is an observable, Eq. ( |l~7| ) becomes a third- 
order polynomial regression formula for [Az(z)] 2 ; 
and one, moreover, with no linear term. 

Fig. 9 shows a plot of [Az(z)] 2 versus (1 + z), 
based on published spectra of 21 different quasars 
with a total of 1320 Lya absorption lines, weighted 
according to resolution (see [^] for details). At first 
sight it does not seem that the regression curve is 
strongly constrained by the data. However, the fit 
is in fact remarkably robust. The reason for this 
is that the curve can consist of only three compo- 
nents: a constant, a downward-opening quadratic, 
and a cubic originating at (0,0). The regression 
coefficients that meet these conditions span only a 
very narrow range of values, and lead directly to 
the 2a (statistical) results |8( 



(18) 



Q M , = 0.014 ±0.006 
ft A .o = 1.08 ± 0.03 . 



similar analysis has been performed under the assumption 
that Lya absorbers make up a homogeneous population of 
clouds with constant comoving size, and it leads to results 
consistent with those presented here ]62|. 
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Fig. 9. Best fit curve of Az(z) to the data from Lya ab- 
sorption lines in the redshift range 1.8 ^ z ^ 4.5. Data are 
weighted according to the quality of the spectra. Dashed lines 
show the quadratic and cubic components of the fit. The co- 
efficients a , a 2 and a 3 are proportional to JIa.o, (1 — SItot.o) 
and Om.o respectively. 



Eqs. ((T^) define what we have referred to as the 
Abar model; this term could equally be applied 
to other vacuum-dominated models without signif- 
icant amounts of CDM (e.g., Fig. 10). 

The Abar model presented in Eqs. (|i~8| ) passes 
several basic consistency tests. Firstly, the sum of 
fi Mj0 + Q Aj0 matches that seen in the CMB exper- 
iments (§|6|). Secondly, the value of f2 Mi0 is within 
the bounds imposed by cosmic nucleosynthesis (§|3|), 
favoring low values of 0, BAR and high values of h 
(1W «S O.OWh- 2 for h = 0.9, or h ^ 0.71 for 
Q BAK = 0.01V 2 )- And thirdly, the regression curve 
(heavy solid line in Fig. 9) passes through z = 
at Az 0.009, in excellent accord with the dis- 
tribution of galaxy structure seen in our own cos- 
mic neighborhood |5(j| (the empty rectangle in the 
upper left-hand corner of Fig. 9). These phenom- 
ena involve independent physics on widely different 
scales, and we would regard it as remarkable for a 
simple procedure like the one described above to 
agree with all three by chance alone. 

The third point however deserves some clarifica- 
tion. The fact that the scale of the local galaxy 
distribution matches that of Lya absorbers agrees 
with our simple picture, in which both types of 
matter cluster predominantly along the cell walls 
(Fig. 8). However, there is still significant debate 
in the literature over the extent of correlations be- 



tween Lya absorbers and other structures [101]. 
Why not "calibrate" the regression curve at z = 
using absorbers rather than visible galaxies? Ultra- 
violet Lya spectra are now available from the Hub- 
ble Space Telescope (Hst) over the redshift range 
^ z ^ 1.3 Q. The absorber material, however, is 
not spread uniformly over the cell walls, but rather 
occurs in filaments and knots. Lines of sight to a 
few quasars, therefore, are far less likely to "no- 
tice" the cell walls at low redshifts, where cell size 
is enormous. Indeed, the Hst spectra show a mean 
redshift spacing < Az >~ 0.04, four times larger 
than that seen in the galaxy distribution [|5C|j . In 
the local universe, in other words, galaxy surveys 
are better guides to large-scale structure than Lya 
absorption spectra. 

A different restriction comes into play at high 
redshifts, where one might expect from the above 
argument that Lya spectra would be increasingly 
reliable. It is certainly true that lines of sight to 
quasars at z > 4 will "notice" most of the cell walls 
they intersect, since the cells themselves are smaller 
(by a factor 1+z) at these redshifts. Another effect, 
however, also goes as 1 + z: Doppler broadening of 
lines due to peculiar motions of the Lya absorbers 
in the cell walls. At z ~ 4 the latter will already be 
of order ~1500 km/s rather than ~300 km/s, intro- 
ducing spurious lines which would masquerade as 
small-scale structure. Recent Keck/HiRES spectra 
of a z = 4.1 quasar confirm this suspicion, yielding 
a mean redshift spacing <Az>~ 0.002 flgj}] , three 
times smaller than the fit to the regression curve 
at z ~ 4 (Fig. 9). The Lya method outlined here 
is most sensitive to the cosmological parameters for 
quasars in the redshift range 2 < z < 4. The most 
recent Keck/HiRES spectra of a quasar at z = 3.1, 
for example, shows <Az>~ 0.006 [fnj], in excellent 
agreement with Fig. 9 at z ~ 3. 

We have not addressed the possibility of bias 
due to inherent evolution within the cell structure. 
There is still debate about the mean comoving scale 
of the Lya distribution in the literature, to say 
nothing of its possible time rate of change. To ob- 
tain higher values of f2 M ,o from the above analysis, 
however, evolution must be in the right direction. 
An appeal to ionization, for instance, should not 
ionize high-redshift clouds more than local ones; 
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this would raise Az at high z and lower the in- 
ferred value of f2 M ,o- Higher matter densities can in 
principle be obtained if evolution is such that sub- 
structure increases with redshift. Observations and 
numerical simulations do suggest the possibility of 
such a trend in the Lya forest, as sheet-like struc- 
tures give way to filaments and knots with time 



[101]. It has been estimated that allowance for an 
effect of this kind might raise the value of fi M ,o horn 

0.05 



120] 



that in Eq. ( J18f) to as much as 

Questions of a more mathematical nature may 
be raised by the smallness of the uncertainties in 
Eqs. (|l^). The possible impact of spectral reso- 
lution, equivalent width, and line blending on the 
line-counting procedure have been considered in 
pC| ], with the conclusion that corrections arising 
from these factors will be minor. The statistical 
robustness of the Lya method can be attributed to 
the absence of a linear term in Eq. (||) for Hubble's 
parameter. More work, however, could be done to 
reduce the possibility of unmodelled systematic er- 
rors. Ultimately one would like to see Az(z) ex- 
tracted from a power spectrum analysis using high 
signal-to-noise spectra (like those now coming from 
Keck/HiRES) together with automated line-fitting 



and counting procedures |7l| • 

We move now to the second argument for a Abar- 
type universe, one which relies on a new analysis of 
the angular power spectrum of the CMB (Fig. 10). 
While the angular positions of the peaks in this 
spectrum fix the sum of matter and vacuum densi- 
ties (§||), their relative heights are largely a function 
of the matter density alone. The odd-numbered 
peaks are produced by regions of the primordial 
plasma which have been maximally compressed by 
infalling material, while the even ones correspond to 
maximally rarefied regions which have rebounded 
due to photon pressure. A high baryon-to-photon 
ratio enhances the compressions and retards the 
rarefractions, thus suppressing the size of, e.g., the 
second peak relative to the first. The strength of 
this effect depends on the fraction of baryons (rela- 
tive to the more weakly-bound neutrinos and CDM 
particles) in the overdense regions. 

Data taken by the Boomerang and Maxima ex- 
periments appear to show an almost total suppres- 
sion of the second peak relative to the first, incon- 




Fig. 10. Observational data on the CMB power spectrum 
from the Boomerang (filled circles) and Maxima experi- 
ments (open circles), together with theoretical expectations 
based on a ACDM model (left) with f2 B AR = 0.039, fl CD M = 
0.317, f2 A ,o = 0.644 and h = 0.7; and a closed Abar model 
(right) with D BA r = 0.034, Ocdm = 0, fi A .o = 1.006 and 
h = 0.75. Figure courtesy S. McGaugh 0. 



sistent (at the 99% level) with expectations based 
on the ACDM model (Fig. 10, left-hand side). The 
ratio of baryons to CDM in the primordial plasma, 
therefore, appears to be higher than predicted. A 
first reaction might be to keep the CDM and raise 
the baryon density; this however brings the theory 
into immediate conflict with nucleosynthesis limits 
on rj BAR |113(] . Other remedies (within the frame- 
work of gravitational instability theory) include tilt- 
ing the spectrum of initial perturbations to disfa- 
vor smaller-scale (higher-order) peaks |3^], erasing 
these peaks outright with processes such as delayed 
recombination [Bq] or decoherence, and varying one 



or more constants of nature [124]. 

The alternative is to take the apparent lack of 
CDM at face value. This can either be done in 
a half-hearted or whole-hearted way. The half- 
hearted way is to retain a minimum density of 
CDM with a statistical "prior." Thus, requiring 
that f^cDM > 0.1, but otherwise fitting the com- 
bined Boomerang, Maxima and Cobe data, one 
obtains a model with best-fit parameters B ar = 
0.032V 2 and CDM = 0.14V 2 @- 

The whole-hearted approach, which may however 
require extending the standard picture of struc- 
ture formation, is to drop the requirement of CDM 
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altogether. Results are shown in Fig. 10 (right- 
hand side), which is a statistical fit to BAR with 
H CDM = and h = 0.75 [35|. The values of Q Afi 
and AT (amplitude) are fixed by the position and 
height of the first peak. The best-fit model passes 
neatly through both peaks (x 2 = 0.85) and has 



n 



M,0 — "BAR — 0.034 

n Ai0 = 1.006 , 



(19) 



in good agreement with primordial nucleosynthe- 
sis (§Ej|) as well as Eqs. (0). This model has an 
age of 22.2 Gyr (with h = 0.75), and a total 
density parameter slightly above one, in agreement 
with suggestions from other new analyses of the 
Boomerang and Maxima data If 



124fl . The 



shape of the CMB power spectrum, along with the 
analysis of Lya absorption lines presented above, 
thus favors an old, closed Abar model, dominated 
by vacuum energy, with no significant contributions 
from CDM or neutrinos — a universe composed al- 
most exclusively of "fire" and "earth." 

8 Conclusions 



both play far more significant roles in determin- 
ing the past and future evolution of the Universe, 
though this is not certain. What is clear is that all 
three forms of matter are dwarfed in importance by 
a newcomer whose physical origin remains shrouded 
in obscurity: the energy of the vacuum. 

We have devoted the last part of our review to the 
hypothesis that vacuum energy dominates so com- 
pletely that there is no room for significant amounts 
of neutrino or exotic dark matter at all. This would 
considerably simplify our picture of the Universe, 
ease problems with the "preposterously" fine-tuned 
values of the observed cosmological parameters, and 
allow more time for galaxies and other structures to 
form. It would also, however, require that we mod- 
ify the GI paradigm. We have reviewed the var- 
ious lines of observational argument, both for and 
against such an idea. It appears to us quite possible 
that the vacuum density f2 Aj0 is close to one, that 
the sole contributions to the matter density £7 Mj o 
come from a small amount of baryons and neutri- 
nos, and that H A>0 and f2 M)0 together are enough to 
"just close" the Universe. 



We have reviewed the evidence for the "four ele- 
ments" of modern cosmology: baryons ("earth"), 
exotic cold dark matter ("water"), neutrinos and 
photons ( "air" ) and vacuum energy ( "fire" ) . Recall- 
ing Eqs. (|), (0), (U), dl|) and (H), we may sum- 
marize the present contributions of each element to 
the total density of the Universe with the following 
ranges of values: 



9-T 



n 



CDM 



TOT.O 



0.012 - 0.041 
" 0.1-0.5 (GI theory) 
- 0.4 (otherwise) 
0.0003 - 0.2ft CDM 
" 0.3-1.1 (GI theory) 
0.4-1.2 (otherwise) 



0.99 - 1.24 [68] , 



(20) 



where "GI" refers to gravitational instability the- 
ory. Baryons, the stuff of which we are made, are 
apparently little more than a cosmic afterthought. 
Neutrinos and the elusive cold dark matter may 
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Appendix: Axions and WIMPs 

This appendix is a brief introduction to axions and 
supersymmetric weakly interacting particles (WIMPs); 
and in particular, to their main claim as compelling cold 
dark matter (CDM) candidates: the theoretical expecta- 
tion that either one, if it exists, naturally has a collective 
density fl C r>M close to the critical density. 

Axions come about as part of the solution to the 
"strong CP problem" in quantum chromodynamics 
(QCD). This theory has had some enormous successes, 
notably in accounting for quark confinement and asymp- 
totic freedom. However, it contains a dimensionless free 
parameter which must be less than ~ 10~ 10 in order 
to explain why strong interactions do not violate par- 
ity (P) or charge-parity (CP). (Upper bounds on the 
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electric dipole moment of the neutron tightly constrain 
any such violations.) To remove this unnaturally small 
number from the theory, particle theorists suppose that 
is driven toward zero by the spontaneous breaking of 
a new symmetry of nature at energy scales f a . The ax- 
ion (a) is the new boson which appears as a result of 
this phase transition, eventually acquiring a rest mass 
m a c 2 oc 1/ f a - Unfortunately neither parameter is con- 
strained a priori by theory (although one hopes that f a 
is less than 10 10 times the QCD phase transition energy, 
or little would be gained from the whole mechanism). 
Observation, however, comes to the rescue. 



Most of the axions with ; 



> 10 2 eV are produced 



thermally, and one can show by solving the Boltzmann 
equation that their combined present mass density would 
be n a w (m a c 2 /130 eV)h~ 2 @. Those in the range 
19 eV < m a c 2 < 32 eV would therefore provide fl a « 
0.3, making them potential CDM candidates. However, 
axions decay into photons on timescales proportional to 
m^ 5 . With masses as large as this, they would decay 
rapidly enough to flood the night sky with ultraviolet 
light. Consistency with the observed intensity of the 
extragalactic background light then leads to the upper 
bound m a c 2 < 3 eV Q. 

Thermal axions with 10 -3 eV < m a c 2 < 3 eV can be 
eliminated on different astrophysical grounds: they cou- 
ple so weakly to ordinary matter that they could stream 
more or less freely out of the cores of red giants and su- 
pernovae, taking energy with them. Observations of the 
neutrino flux from supernova SN1987a show no evidence 
for such an effect |77j] . 

Axions with m a c 2 < 10~ 3 eV, finally, are largely 
nonthermal, with a collective density given by f2 a « 



(m a c 2 /4:ficV)- 7 ^h- 2 |pl0|.Fl Rest masses 



c 2 < 

lOfieV are then disqualified because they would provide 
too much CDM (Cl a > 0.5). The axion is thereby re- 
stricted to a relatively small window of potential rest 
masses, 10 fieV < m a c 2 < 1000 fxeV. Its plausibility as 
a CDM candidate rests on the fact that this allowed 
window encompasses the range of values (15 fieV < 
m aC 2 < 25 /xeV) corresponding to Q a w 0.3. The slow 
decay rate of axions in this mass range can, moreover, 
be enhanced by trapping them in strong magnetic fields; 
this is the basis for a number of ongoing experimental 
detection efforts in Japan and the U.S.A. [ 1 10 1 . 

Like axions, supersymmetric WIMPs have their ori- 
gin in a new symmetry of nature, spontaneously broken 
in the early Universe. This is supersymmetry (SUSY), 
which pairs each boson with a fermion superpartner, and 
vice versa. Because no such pairs can be formed with 
the known bosons and fermions of the standard model 



of particle physics, the number of fundamental particles 
must be doubled. The new SUSY partners are presum- 
ably so massive that they have not been discovered yet; 
the lightest SUSY particle (LSP) in particular must have 
a rest mass m x c 2 > 50 GeV |36| . 

This LSP plays the role of the WIMP in SUSY the- 
ories. It is stable, decaying solely by pair-annihilation 
with itself (a very slow process). This is due to an ad- 
ditional new symmetry of nature, known as i?-parity, 
which is necessary (in "minimal SUSY" models) to keep 
the proton from decaying via intermediate SUSY states. 
(There are also nonminimal SUSY theories in which this 
symmetry too is spontaneously broken, and the proton 
can decay.) i?-parity requires the number of SUSY (and 
non-SUSY) partners to be conserved in any reaction, so 
that, as the Universe cools, heavier SUSY particles can 
break down into lighter ones, but not into ordinary parti- 
cles. Eventually, most of the SUSY mass in the Universe 
thus ends up in the form of LSPs. 

Using the Boltzmann equation, one can calculate the 
collective density fl x of these particles in terms of a num- 
ber of free parameters such as the SUSY-breaking energy 
scale, and the composition of the LSP.0 Because its 
annihilation cross-section goes as w~ 2 , one finds that 
an LSP more massive than ~ 3 TeV would annihi- 
late so slowly as to overdose the Universe. This leaves 
an available SUSY WIMP mass window of 50 GeV < 
m x c 2 < 3 TeV. The collective LSP density turns out to 
lie within three orders of magnitude of the critical den- 
sity over most of this range [ 1 17 1 , which makes the SUSY 
WIMP, like the axion, a plausible CDM candidate. 

Many experiments around the world are currently 
searching for WIMPs in this mass range, assuming for 
example that they are gravitationally bound in the 
galactic halo and will occasionally be responsible for 
scattering events in target nuclei as the Earth follows 
the Sun around the Milky Way O] . One team (Dama) 



has claimed evidence for a 



= 59+11 GeV WIMP 



signature at 3tr confidence |T3| , but this is disputed by a 
second (Cdms) which has searched a larger region of pa- 
rameter space and found nothing Q| . Strong constraints 
have also been placed on WIMPs in nonminimal versions 
of SUSY, where the LSP can decay into photons and con- 
tribute excessively to the intensity of the diffuse x- and 
7-ray backgrounds p2] |. We watch these developments 
with interest. 



1 This is under debate, and may be up to an order of mag- 
nitude larger if string effects play an important role Jilt . 



19 This is most likely to be the neutralino, a linear superpo- 
sition the photino, the zino and two higgsinos (the spin-1/2 
SUSY partners of the photon, Z and Higgs bosons). A less 
favored candidate, because it annihilates too slowly, is the 
gravitino, the spin-3/2 SUSY partner of the graviton ]35[|. 
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